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Abstract 


In  this  dissertation,  two  nondestructive  dual  ridged  waveguide  (DRWG)  material 
characterization  systems  are  investigated.  The  single  and  clamped  DRWG  probe 
geometries  were  analyzed  in  previous  work;  however,  that  research  only  incorporated 
the  dominant  DRWG  mode.  Here,  that  restriction  is  removed  and  the  existence  of 
evanescent  higher-order  modes  is  considered.  Theoretical  analysis  of  the  single  and 
clamped  DRWG  probes  is  presented  and  discussed.  The  approach  taken  here  is  similar 
to  that  presented  in  previous  research:  Love’s  equivalence  theorem  is  used  to  remove  the 
DRWG  apertures  which  are  replaced  with  equivalent  magnetic  currents  radiating  in  the 
presence  of  the  background  parallel-plate  waveguide  structure.  Enforcing  the  continuity  of 
the  tangential  magnetic  fields  in  the  DRWG  and  parallel-plate  regions  yields  a  system  of 
coupled  magnetic  field  integral  equations  (MEIEs).  This  coupled  MEIE  system  is  solved 
using  the  Method  of  Moments  (MoM)  where  the  tangential  electric  and  magnetic  fields 
in  the  DRWG  are  used  as  expansion  and  testing  functions,  respectively.  Inversion  of  the 
resulting  MoM  impedance  matrix  produces  theoretical  expressions  for  the  reflection  and 
transmission  coefficients  which  are  numerically  inverted  (via  nonlinear  least  squares)  to 
yield  estimates  of  the  permittivity  and  permeability  of  the  material  under  test.  While 
the  steps  taken  in  the  theoretical  development  are  similar  to  previous  work,  the  addition 
of  higher-order  modes  into  the  analysis  substantially  complicates  the  derivation  and  is 
a  significant  extension  of  the  existing  dominant-mode-only  literature.  Lastly,  simulation 
results  of  the  two  structures  are  presented.  A  comparison  of  the  dominant-mode  only 
reflection  and  transmission  coefficients  with  the  higher-order  mode  coefficients  is  provided. 
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MULTI-MODE  ANALYSIS  OF 


DUAL  RIDGED  WAVEGUIDE  SYSTEMS 
FOR  MATERIAL  CHARACTERIZATION 


I.  Introduction 


1.1  Introduction 

Electromagnetic  material  characterization  is  the  process  of  determining  the  dielectric, 
magnetic,  and  magnetoelectric  properties  of  a  material.  For  simple  (i.e.,  linear, 
homogeneous  and  isotropic)  media,  which  is  the  focus  of  this  research,  these  properties 
are  the  permittivity  e  =  ere  +  j^im  and  permeability  jj.  =  Hre  +  In  general,  material 
characterization  is  classified  as  destructive  or  nondestructive.  In  destructive  techniques, 
the  material  under  test  (MUT)  is  machined  or  altered  in  some  way  in  order  to  fit  into  the 
electromagnetic  testing  device  (e.g.,  waveguide).  Destructive  techniques  often  have  the 
significant  advantage  of  allowing  extraction  of  the  permittivity  and  permeability  of  a  MUT 
via  simple  closed-form  analytic  expressions  which  are  computationally  efficient  (e.g.,  1-2 
seconds).  However,  since  destructive  techniques  require  the  MUT  to  be  cut  to  fit  inside  the 
test  device,  significant  errors  can  be  introduced  if  sample  fit  is  not  perfect.  In  contrast,  a 
nondestructive  process  seeks  to  determine  material  parameters  without  altering  the  shape 
of  the  MUT.  Thus,  sample  fit  problems  are  mitigated  and  substantial  time  is  saved  in  not 
having  to  perform  any  sample  preparation.  In  addition,  nondestructive  techniques  have  the 
benefit  to  be  utilized  in  the  field  as  point  inspection  tools.  The  drawback  of  nondestructive 
techniques  is  that  the  measurement  environment  is  substantially  more  complex,  thus 
computationally-intensive  iterative  algorithms  are  required  to  extract  permittivity  and 
permeability  values  of  the  MUT.  Due  to  the  intended  purpose  of  developing  fiight-line 
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testing  tools,  the  focus  of  this  research  will  be  on  nondestructive  evaluation  of  simple 
media.  While  many  choices  for  nondestructive  evaluation  (NDE)  exist,  this  research  will 
focus  and  expand  upon  prior  research  using  Dual  Ridged  Waveguide  (DRWG)  probes. 

1.2  Problem  Statement 

NDE  of  material  parameters  has  always  been  important  to  the  United  States  Air 
Eorce  (USAE).  Material  coatings  are  predominant  on  many  air  vehicles  to  manage  radar 
detectability,  thus  material  characterization  on  the  flight  line  plays  a  crucial  role  in  vehicle 
health  performance  and  maintenance.  Since  these  material  coatings  are  physically  adhered 
to  the  vehicle  body  (i.e.,  metallic  conductor),  nondestructive  techniques  are  required  to 
assess  permittivity  and  permeability  data,  and  thus,  the  radar  reducing  performance  metrics. 
In  order  to  accurately  locate  and  detect  permittivity  and  permeability  changes  in  the  aircraft 
material  coating,  NDE  test  fixtures  having  large  bandwidth  are  desirable  due  to  enhanced 
resolution.  Current  techniques,  namely  single  and  double  rectangular  waveguide  probes 
used  to  characterize  perfect  electrical  conductor  (PEC)  backed  material,  lack  appropriate 
bandwidth.  A  clamped  rectangular  waveguide  technique  useful  in  laboratory  environment 
also  exists,  but  cannot  be  utilized  for  characterizing  conductor-backed  materials  and  suffers 
the  same  bandwidth  limitation. 

The  primary  goals  of  this  research  effort  are  to  replace  the  existing  rectangular 
waveguide  probes  (single,  double  and  clamped)  with  dual-ridged  waveguide  components. 
The  dramatic  benefit  of  using  dual-ridged  waveguides  is  a  three-fold  improvement  in 
bandwidth  (and  hence,  three-fold  improvement  in  resolution).  This  effort  will  build 
upon  other  work  related  to  dual-ridged  waveguides,  namely,  single  and  clamped  probe 
research  performed  by  [19,  20],  which  considered  dominant-mode  propagation  only. 
The  investigation  here  considers  the  existence  of  higher-order  modes  in  the  dual-ridged 
waveguides  and  how  they  influence  the  accuracy  of  extracted  permittivity  and  permeability 
values.  This  is  a  critical  issue  since  the  inclusion  of  higher-order  modes  can  add  substantial 
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time  to  the  material  property  extraction  process,  namely,  from  several  minutes  for  the 
dominant  mode  to  several  days  for  the  inclusion  of  several  higher-order  modes.  Obviously, 
a  material  extraction  process  that  takes  several  days  is  quite  useless  for  the  high  demands 
of  a  flight  line.  If  it  can  be  shown  that  the  higher-order  modes  do  not  significantly 
contribute  to  accuracy,  than  the  dual-ridged  waveguide  can  be  considered  a  valuable  and 
time-efficient  diagnostic  tool,  at  least  for  the  class  of  materials  typically  encountered  in 
night-line  environments.  In  addition  to  the  use  of  DRWG’s  for  single  and  clamped  probes, 
a  double  (i.e.,  dual)  probe  geometry  will  also  be  investigated.  The  dual  probe  dual-ridged 
waveguide  test  fixture  has  the  advantage  of  being  able  to  acquire  both  reflection  and 
transmission  data,  which  is  important  in  material  characterization  since  these  quantities 
are  linearly  independent. 

Investigation  of  literature  reveals  that  the  development  of  the  dual  probe  DRWG 
system  as  well  as  higher  order  mode  expansion  of  the  single  and  clamped  dual-ridged 
systems  represent  original  work.  They  also  have  the  potential  of  being  considered  as 
significant  contributions  as  mentioned  above,  especially  for  flight-line  applications.  This 
research  may  also  lead  to  better,  more  accurate,  NDE  techniques  for  broader  classes  of 
media. 

1.3  Material  Characterization  Background  and  Methods 

It  is  important  to  understand  the  benefits  of  accurate  NDE  of  materials.  As 
advancements  in  materials  and  material  manufacturing  are  made,  accurate  characterization 
is  important  for  the  intended  application.  As  one  colleague  wrote,  “when  producing  a 
material,  how  do  you  know  you  made  what  you  wanted  to  make?”  [30].  Without  being 
able  to  confirm  a  material’s  properties,  it  would  not  be  possible  to  accurately  model  new 
designs  such  as  antennas  and  absorbing  coatings.  While  many  numerous  techniques  exist  in 
the  microwave  spectrum,  their  effectiveness  varies.  Material  type,  models  used,  frequency 
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range,  and  lab  set-up  all  play  a  role  in  determining  whieh  teehnique  is  most  effeetive.  This 
section  presents  a  short  review  of  several  techniques  involving  waveguides. 

1.3.1  Coaxial  Waveguides. 

Coaxial  waveguides  are  very  common  due  to  their  bandwidth  and  information¬ 
carrying  capabilities.  These  waveguides  consist  of  an  inner  conductor  surrounded  by 
insulation  which  is  further  surrounded  by  an  outer  conductor  sometimes  referred  to  as  a 
shield.  This  is  further  surrounded  by  another  layer  of  protective  insulation.  A  picture  of  a 
generic  coaxial  cable  structure  is  shown  in  Figure  (1.1).  These  cables  are  generally  round 
and  flexible  and  are  commonly  used  for  material  characterization  when  the  region  between 
the  inner  and  outer  conductors  is  air  filled  (center  conductor  held  into  position  by  small/thin 
teflon  dielectric  beads). 


COPPER 

WIRE 


R  MESH 
UTSIDE  INSULATION 


Figure  1.1:  Coaxial  cable  structure. 


Coaxial  probes  have  been  used  in  many  applications  [2-4,  6,  8,  11,  12,  17,  25,  29,  32, 
36].  While  they  provide  the  potential  for  broadband  measurements,  these  systems  depend 
greatly  on  the  MUT  and  waveguide  size  in  order  to  obtain  accurate  material  properties 
[11,  12].  Many  of  these  references  discuss  that,  although  the  theoretical  bandwidth  of 
a  commonly-used  7mm  coax  is  18  GHz,  the  effective  bandwidth  for  accurate  material 
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characterization  is  only  about  10  GHz  due  to  the  field  strueture  and  laek  of  bore-site  field 
eoupling  into  the  material  under  test,  espeeially  at  lower  frequeneies.  It  has  also  been  shown 
that  field  penetration  is  better  in  lossy  material  when  eompared  to  lowloss  materials,  again 
due  to  field  strueture  and  the  types  of  eurrents  that  ean  exist  in  these  types  of  materials. 

1.3.2  Rectangular  and  Circular  Waveguides. 

Reetangular  and  eireular  waveguides  are  mueh  more  physieally  robust  than  eoaxial 
waveguides.  There  eonstruetion  is  also  mueh  simpler  being  a  simple  reetangular  or  eireular 
tube  eonstrueted  of  a  material  with  near  PEC  properties.  A  reetangular  and  eireular 
waveguide  are  shown  in  Figure  (1.2). 


Figure  1.2:  Reetangular  and  eireular  waveguide. 
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Rectangular  and  circular  waveguides  have  also  been  used  in  many  applications 
[7,  9,  18,  19,  23,  24,  26,  31,  34].  In  addition  to  there  robustness,  rectangular  and  circular 
waveguides  do  not  have  the  coupling  problems  previously  noted  for  coaxial  waveguides. 
They  also  have  very  simple  field  expressions  in  terms  of  elementary  functions  (sines, 
cosines,  exponentials  and  Bessel  functions)  and  corresponding  propagation  constants  of  the 
allowed  modes  can  be  found  in  closed  form.  The  primary  drawback  of  these  test  fixtures  is 
that  they  are  band  limited.  As  an  example,  an  air-filled  X-band  waveguide  (having  cross- 
sectional  dimensions  of  2.286  cm  by  1.016  cm)  has  a  bandwidth  of  only  4  GHz  when 
operated  in  the  traditional  dominant-mode  configuration. 

1.3.3  DRWGs. 

DRWGs  are  somewhat  similar  to  rectangular  waveguides  in  construction.  They 
resemble  a  rectangular  waveguide  with  a  ridge  on  the  long  sides  (typically  considered 
the  X  axis)  of  the  waveguide.  This  gives  the  waveguide  the  shape  of  an  H.  Figure  2.1 
shows  a  cross-sectional  depiction  of  a  dual  ridge  waveguide  (single  ridge  waveguides  are 
also  common,  but  are  not  investigated  here  due  to  their  limited  bandwidth).  Originally 
developed  in  [27],  these  waveguides  have  the  benefit  of  robustness  of  a  solid  waveguide 
and  a  three-fold  increase  in  bandwidth  over  a  similar  rectangular  waveguide.  However,  the 
field  expressions  are  more  complicated  as  a  result  of  the  three  distinct  regions  (i.e.,  separate 
field  expressions  for  each  region).  In  addition,  the  allowed  propagation  constants  must  be 
found  numerically,  leading  to  reduced  computational  efficiency.  Recent  research  efforts 
have  been  conducted  on  single  probe  DRWG  systems  [21]  as  well  as  clamped  systems  [20] 
showing  the  ability  to  accurately  characterize  simple  materials  with  approximately  3  times 
the  bandwidth  of  a  traditional  X-band  waveguide. 

1.3.4  Single  Probe  Waveguide  Systems. 

A  single  probe  system  uses  only  a  single  waveguide  to  characterize  a  material. 
The  waveguide  can  be  of  any  type  previously  discussed.  A  single  probe  waveguide 
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system  is  one  of  the  subjects  for  this  research  and  an  example  can  be  seen  in  Chapter 
3.  While  numerous  papers  have  been  published  on  single  probe  waveguide  systems 
[7-9,  11,  12,  17,  18,  25,  26,  29,  31,  32,  34,  36],  a  single  probe  system  cannot  easily 
determine  both  complex  permittivity  and  permeability.  This  is  because  the  system  can 
only  take  a  single  measurement  (i.e.,  reflection  measurement)  allowing  for  calculation  of  a 
single  unknown.  Several  techniques  have  been  developed  to  alleviate  this  problem  such 
as  two  thickness  method  [9,  10,  26,  33],  two  layer  method  [3,  14,  26,  34],  frequency 
varying  method  [26,  35],  free  space  backed  method  [3,  25,  34],  and  two  iris  method 
[13].  However,  these  techniques  have  disadvantages  such  as  instability  (two  layer  method), 
physically  difficult  to  measure  (free  space  backed  method),  require  prior  knowledge  of  the 
MUT  (two  thickness  and  frequency  varying  methods),  or  require  two  waveguides  (two 
iris  method)  [20].  As  mentioned  previously,  the  accommodation  of  higher-order  modes 
in  a  dual-ridged  single  probe  waveguide  for  the  purpose  of  enhancing  material  property 
extraction  is  investigated  as  part  of  this  research  effort,  and  is  discussed  in  Chapter  3. 

1.3.5  Two  Probe  Waveguide  Systems. 

A  two  probe  system  uses  two  separate  waveguides  typically  in  a  clamped  [2,  24] 
or  dual  [6,  23]  configuration.  The  waveguide  can  be  of  any  type  previously  discussed. 
Dual-ridged  waveguide  dual  probe  and  clamped  systems  are  the  remaining  effort  of 
this  research  and  are  discussed  further  in  Chapters  2  and  4,  respectively.  As  noted 
previously,  two  probe  systems  are  desirable  because  they  can  make  multiple  independent 
measurements  (typically  reflection  and  transmission)  and  thus  determine  both  permittivity 
and  permeability  simultaneously. 

1.4  Assumptions 

The  following  established  assumptions  (e.g.,  [30])  are  used  in  this  research. 

•  An  time  dependence  is  assumed  and  suppressed  throughout. 
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•  The  transverse  axes  are  x  and  y  and  the  direction  of  propagation  is  the  z  axis. 

•  Conducting  surfaces,  such  as  the  waveguide  walls  and  parallel  plates,  are  assumed  to 
be  a  PEC. 

•  The  transverse  dimensions  of  the  material  sample  and  parallel  plates  are  assumed  to 
be  infinite  in  extent. 

•  The  DRWG’s  contain  only  free  space  with  material  parameters  6o  and  //q- 

•  The  MUT  is  simple  (linear,  isotropic,  homogeneous),  and  the  sample  has  a  uniform 
thickness. 

•  The  waveguide  probes  are  perfectly  aligned  in  all  dimensions. 

1.5  Overview  and  Organization 

This  chapter  has  provided  the  background  and  motivation  for  this  research.  Several 
NDE  techniques  were  briefly  discussed  along  with  previous  work  on  the  subject.  Goals  for 
this  research  as  well  as  how  this  research  will  impact  the  USAE  was  presented.  Chapter 
2  will  discuss  the  dual  probe  DRWG  system  development.  Chapters  3  and  4  will  cover 
higher  order  expansion  of  the  single  probe  system  and  clamped  probe  systems,  respectively. 
Einally,  Chapter  5  will  present  conclusions  of  this  work  and  suggest  future  research  areas. 
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11.  Dual  Probe  DRWG 


In  this  chapter,  the  analysis  of  a  dual  probe  dual-ridged  waveguide  (DRWG),  shown 
in  Figures  2.1  and  2.2,  is  performed.  The  analysis  includes  an  examination  of  the 
reflection  and  transmission  (i.e.,  coupling)  coefficients  of  the  dual  probe  DRWG  placed 
on  a  conductor-backed  material.  It  is  assumed  in  the  analysis  that  only  a  single  dominant 
mode  is  propagating  in  the  free-space  filled  DRWG’s.  Two  simulation  scenarios  are 
considered,  namely,  a  lossy  and  low  loss  conductor-backed  material.  The  fundamental 
goal  for  performing  this  analysis  is  to  assess  if  this  probe  has  sufficient  sensitivity  to 
simultaneously  extract  permittivity  and  permeability  of  the  material,  and  thus,  comprise 
a  valuable  flight-line  diagnostic  tool.  This  can  be  answered  through  the  investigation  of  the 
reflection  and  transmission  coefficient  amplitudes.  Recall,  if  two  material  properties  are 
desired,  then  two  linearly-independent  measurements  are  required  (in  this  case,  reflection 
and  transmission  coefficients).  However,  it  is  also  necessary  that  the  measured  coefficients 
have  sufficient  signal-to-noise  ratio.  The  dual  probe  analysis  aims  to  answer  this  critical 
question.  It  is  intuitively  expected  for  the  lossy  case  that  transmission  should  be  maximized 
(which  is  desired)  when  dual  probe  separation  is  minimized  due  to  energy  absorption  in  the 
conductor-backed  material  as  well  as  wave  spreading.  However,  will  minimized  probe 
spacing  be  sufficient  for  traditional  radar  absorbing  materials?  The  intended  analysis  will 
answer  this  critical  question.  In  regards  to  lossless  or  low  loss  media,  is  minimized  probe 
spacing  necessarily  desired?  The  analysis  that  follows  will  demonstrate  that  this  is  not 
necessarily  true  for  this  case. 

The  dual  probe  DRWG  analysis  will  begin  with  a  brief  review  of  the  dual-ridged 
waveguide  field  expressions.  Field  plots  of  the  dominant  mode  is  provided,  to  demonstrate 
why  sufficient  transmission  coupling  is  a  concern.  Next,  the  fields  in  the  parallel-plate 
region  of  the  dual  probe  DRWG  are  found  with  the  aid  of  Love’s  Equivalence  Theorem. 
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Subsequent  enforeement  of  tangential  boundary  conditions  at  the  dual  probe  DRWG 
apertures  leads  to  coupled  magnetic  field  integral  equations  (MFIE’s)  for  the  unknown 
reflection  and  transmission  coefficients.  These  MFIE’s  are  solved  via  the  Method-of- 
Moments  for  both  a  lossy  and  low  loss  material.  Details  of  the  complex-plane  analysis 
utilized  in  the  solution  process  are  included  for  the  benefit  of  the  reader.  Simulation 
results  for  the  reflection  and  transmission  coefficients  of  both  materials  are  presented  and 
discussed,  and  general  conclusions  are  provided. 


Figure  2.1:  Cross-sectional  view  of  dual  probe  DRWG  geometry 
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2.1  DRWG  Fields 

Figure  2.1  shows  the  cross-sectional  view  of  the  free-space  filled  dual  probe  DRWG, 
with  each  DRWG  comprised  of  a  right  trough  region  (A.r<v:<a/2,-b/2<3;<b/2),  a  central  gap 
region  (-Ax<x<Ax,-^y<y<^y),  and  left  trough  region  (-a/2<x<-Ax,-i>/2<y<b/2).  Figure 

2.2  depicts  the  side  view  of  the  dual  probe  DRWG  in  contact  with  the  conductor-backed 
material  having  permittivity  e,  permeability  ju  and  occupying  the  parallel-plate  region  z=0 
to  z=d.  It  is  also  assumed  that  the  adjacent  DRWG’s  are  located  a  distance  of  yo  apart  from 
center  to  center. 

The  following  desired  fields  for  a  DRWG  were  derived  in  [15,  27].  These  fields  are 
required  and  used  throughout  this  research  as  the  fields  in  the  DRWG  regions  of  all  the 
DRWG  systems. 
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The  T transverse  DRWG  field  expressions  are  given  by 

Ex  =  Yj  cos(fcf„;c)  +  (fin  +  an)  sinm„x)]  sin  -  Ay)) 

n 

=  -  y  [{fin  -  an)  sin(fc«„v:)  -{fin  +  an)  eos(fcf„a:)]  eos  (fc^„(y  -  Ay)) 

(2.1) 


■g  f  =  xEl  +  y£f 


= 


zx  e; 


in  the  |jc|  <  Av:  (gap)  subregion, 

K'  =  J]  ane^,n,mCOS 


n,m 


=  Yj  ^’^4n,mSin 


f,,  ,  b\ 

'ii-i 

sin| 

'is  ^  \ 

(2-dl 

eos  1 

kym(y  -  2)] 

in  the  x  >  Ax  (right)  trough  subregion,  and 


y = - 

n,m 

Ey  = 


kl. 


=  xE^^  +  y£; 


sinj^Jy-  -) 


n,m 


E  -  -  X 
^xm  \  2 


eosl^Jy-  -) 


-t‘:  =  xE‘‘  +  yE‘; 


h 


zx  e‘ 
ZTE 


(2.2) 


(2.3) 


in  the  x  <  -Ax  (left)  trough  subregion.  In  the  above  field  expressions,  Z^^  =  (o^iolE  is  the 
T E^  wave  impedanee  and 


f^yn 


eL  = 


2ExnCOS{klnAx) 

1 


2cOS{kinAx) 

kyni4^m,n 


(2.4) 


^  xn,m 


k{n^l{l  +  6„,o)sin[k{J^  -  Av:)] 


xmy  2 


= 

yn,m  ^ 


if/„ 


^(1  +  Smo)sin[k{^{^  -  Av:)] 
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where 


U  1 


^\Kn-k‘y^  J4n  +  k‘y^ 


sin  k‘  -  -  AyU-i-lfsin  kl 


Variables  k\,^  =  mnib  and  kyn  =  rnijlAy  are  the  y-directed  DRWG  wavenumbers  in 
the  trough  and  gap  subregions,  respectively;  k‘^^  and  klm  are  the  a:-directed  DRWG 
wavenumbers  in  the  trough  and  gap  subregions,  respectively;  and  an  and  are  the 
complex  TE^  modal  amplitudes.  Note  that  k^  =  Jkl  -  k^  =  Jkl  -  (k‘  )  -  (k‘^^)  = 


kl  -  {k^yn)  -  {kin)  where  ko  =  a>  ^[e^  and  co  is  the  radian  frequency. 


The  modal  plots  of  the  transverse  electric  field  distribution  for  the  dominant  mode 
of  the  DRWG  is  shown  in  Figure  2.3.  It  is  critical  to  note  that  the  dominant  mode  is 
tightly  confined  to  the  central  gap  region.  This  is  of  great  concern  because  the  effective 
spacing  between  the  adjacent  DRWG’s  energy  concentration  in  the  dual  probe  system  is 
substantially  farther  apart  than  in  a  rectangular  waveguide  dual  probe  configuration.  Thus, 
the  main  question  becomes  is  there  sufficient  transmission  coupling  between  the  DRWG’s 
for  both  a  lossy  and  low  loss  material.  The  goal  of  this  chapter  is  to  answer  that  critical 


concern. 


2.2  Dual  Probe  DRWG  Region  Fields 

As  a  wave  travels  down  the  first  (centered  about  the  origin)  DRWG  and  interacts  with 
the  boundary  between  the  waveguide  and  the  MUT,  a  portion  of  the  energy  will  be  reflected 
back  into  the  waveguide  [5] .  A  portion  of  the  energy  will  also  be  transmitted  into  waveguide 
2  (offset  yo  about  the  origin).  Thus  we  must  first  find  the  total  fields  in  the  DRWG  regions. 
Total  fields  in  the  DRWG  portion  of  the  system  are  found  using  the  DRWG  fields  reviewed 
in  section  2.1.  Note  however,  that  because  of  the  offset  in  the  y  direction  a  shift  is  applied 
to  the  fields  in  the  second  waveguide.  Thus  the  y  variation  of  each  field  component  has 
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-3  for  first  TE^  mode 

X 10  ^ 


mm 


Figure  2.3:  Dominant  mode  transverse  eleetric  field  distribution. 


a  form  similar  to  sin\klm(y  -  (jo  +  A}^))]  where  the  shift  is  given  by  yo-  Allowing  for  this 
yields  fields  in  the  form  of 


(2.6) 

^tDRWG  -  ^  ^  f  (p) 

(2.7) 

^^DRWG  — 

(2.8) 

^tDRWG  -  -^h^Cp) 

(2.9) 

where  F  and  T  are  the  unknown  (i.e.,  to  be  determined)  refleetion  and  transmission 
eoelfioients,  respeetively.  Note,  superscripts  1  and  2  denote  waveguide  1  and  2. 
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2.3  Parallel  Plate  Region  Fields 


As  mentioned  previously,  a  wave  traveling  down  the  DRWG  interaets  with  the 
interfaee  between  the  waveguide  and  the  MUT,  eausing  a  portion  of  the  energy  to  be 
transmitted  into  the  MUT  (i.e.,  parallel  plate  field  region).  This  energy  gives  rise  to 
eleetromagnetie  fields  inside  the  PP  region  of  the  system.  To  determine  these  fields, 
first  Love’s  equivalenee  prineiple  [5]  is  used  to  replaee  the  DRWG  aperture  fields  with 
an  equivalent  souree  in  the  form  of  magnetic  surface  currents,  namely  [5] 

Mi=-fix'Ei.  (2.10) 

M2  =  -hx^2-  (2.11) 

where  h  is  the  normal  unit  vector  pointing  into  the  MUT  region,  with  E  i  and  F  2  to  be 
determined.  This  equivalent  system  is  shown  in  Figure  (2.4).  These  aperture  currents  are 
located  at  z  =  0  and  occupy  the  aperture  field  regions  as  discussed  in  section  2.1. 


Z=0  Z=d 


S,jLl 

\ 

PEC^ 

\PEC 

\ 

\ 

\ 

M. 

\ 

Figure  2.4:  Equivalent  system  for  parallel  plate  region  of  the  dual  probe  DRWG  system 
after  Love’s  Equivalence  Theorem  is  applied. 
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Now  that  the  equivalent  currents  have  been  identified,  a  Vector  Potential  (VP)  solution 
is  used  to  find  the  parallel  plate  region  magnetic  field.  This  solution  is  as  follows  [5]: 


=  —(ei ,  +  V,V-  )FrG^,z).  (2.12) 

jojeix 

where  1_  t  =  xx  +  yy,  and  =  o/en.  However,  the  F  VP  must  first  be  identified.  The 
surface  currents  Mi  and  M2  support  this  VP  through  the  following  relation. 


Ff^(p,z)  =  JJ"  G(p,z[p',0)-eMi(p')ds' 
+  ff  G(p,z[p',0)-eM2(p')ds' 


(2.13) 


where  G(7? ,  z\p  '  ,G)  is  the  parallel  plate  Green’s  function  to  be  discussed  in  the  next  section, 
and  Ml  and  M2  are  found  using  equations  2.10  and  2.1 1,  respectively.  Note,  S 1  and  S2  are 
the  DRWG  aperture  regions  as  previously  discussed.  Equating  the  tangential  magnetic  field 
at  the  apertures  leads  to  the  coupled  MFIE’s 


—{ei  1  +  v,v-  )FrG^,z)  =  (1  -  r)tj(p).  (2.14) 

ja)6n 

t  +  ViV-  )Ff  (p,z)  =  -Tt^Oo)  (2.15) 

ja)6ix 

These  MFIEs  are  solved  in  the  next  section  via  the  MoM. 

2.4  Solving  the  MFIEs 

This  section  details  the  solution  of  the  MFIE’s  of  the  previous  section.  The  first  step 
of  the  MoM  is  to  expand  the  unknown  currents  (Mi)  and  (M2)  in  a  series  of  basis  functions. 
In  this  case  we  will  exploit  the  fact  that  we  know  the  transverse  electric  fields  in  the  DRWG 
region  at  z  =  0“,  and  we  know  that  these  transverse  fields  must  be  continuous  across  the 
boundaries  at  z  =  0"^,  thus  we  can  approximate  the  equivalent  magnetic  currents  as 


Ml  (7^  )—  h  X  E I  ~  h  X  E 


(2.16) 
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M2(p')  =  -hxE2^-hx  (2.17) 

Substitution  of  2.16  and  2.17  into  2.13  and  further  insertion  into  2.14  and  2.15  yields  the 
following  relations: 


+  ff  G(p,z[p',0)-Mi(p')ds'  +  ff  G(p,z[p',0)-M2(p')ds' 

LJJs,  JJS2  J 

=  (i-r)XV) 

(2.18) 

i  +  ViV-)  ff  G(p,zfp',OyMi(p')ds'  +  ff  G(p,zfp',0yM2(p')ds' 

J<^d  [JJs,  JJs2  J 

=  -t"^?(p) 

(2.19) 

Next,  both  sides  of  equations  2.18  and  2.19  are  tested  with  the  eomplex  eonjugate  of 
the  magnetie  field.  The  testing  operators  are  as  follows: 


h}*dSi 


(2.20) 


{•••}•  hfdS: 


for  2.18  and  2.19,  respectively.  After  testing,  2.18  and  2.19  become 


(2.21) 


ff  +  V(V-)  ff  G(p,z[p',0)-Mi(p')ds' 

+  JJ"  G(p,z[p',0)- M2(p')ds'  j-  h}*dS  1 

=  f£  |(i-r)tjG^)j-ti*J5i 

ff  +  V(V-)  ff  G(p,z[p',0)-Mi(p')ds' 

+  JJ"  G(p,z[p',0)- M2(p')ds'  j-  hf*dS2 

=  -J£  lThf(^)jl^f*dS2 


(2.22) 


(2.23) 
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It  is  desirable  to  use  a  speetral  representation  of  the  Green’s  function.  This 
representation  allows  easy  evaluation  of  the  spacial  integrals  involving  v:,  x',y,y',  resulting 
in  simple  trigonometric  functions  as  shown  later.  However,  it  does  add  two  more  integral 
dimensions  over  spectral  variables  r]  and  The  spectral  representation  of  the  Green’s 
function  is  shown  below 

G(^,z[p',0)  =  ^ff_  (2.24) 

where 

~A  =  x^  +  yrj  (2.25) 


and 


with 


G(A,zlO)  =  I 


cosh(pd) 

‘psinh(pd) 


(2.26) 


P  =  if 


(2.27) 


After  transforming  the  Green’s  function,  rearranging  the  order  of  integration,  and 
pulling  out  constants,  equations  2.22  and  2.23  become 

77b  f£  fl,  £  fl  ^  •  (-« 

+  ({fl_  ,  +  V,V  • )  5(1,210)  •  (-h  x'e^^drwg]) 

■  h]*  ds  ds  drj  d^ 

=  ff  ((i-r)t^G^))-ti*  jSi. 

(2.28) 
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‘  rr  rr  rr 

jconiln)^  JJ_oo  JJs2  I  JJs, 

^ffs. 


e  ( {1^1  f  +  V,V  • )  GU, z|0)  •  ( -n  x 


e-j^^  Q-my' 


tDRWG 


{kH_ ,  +  V,V  • )  G(l,  ZIO)  •  ( -n  X 


h^*  ds  ds  dr]  d^ 

h^*dS2 


(2.29) 


Once  all  operations  are  eondueted  on  2.28  and  2.29,  a  series  of  six  dimensional 
integrals  remain  on  the  left  of  the  equalities,  and  a  double  integral  remains  on  the 
right  of  the  equalities.  The  double  integral  to  the  right  has  three  seetions  that  must 
be  ealeulated  separately,  one  for  eaeh  subregion  of  the  waveguide.  These  integrals  are 
straightforward  and  easily  ealeulated  using  standard  trigonometrie  identities  and  integration 
of  trigonometrie  funetions.  The  left  side  integrals  are  more  eomplieated.  In  the  dual 
probe  system  64  total  6  dimensional  integrals  must  be  ealeulated.  One  example  of  these  6 
dimensional  integrals  is  shown  below 


ZZ 

n,m  h,m 


t 

P  p  ~  - 
^yn,m  ^yn,m 


Xoo 

CO  %J  —cx 


cosh(pd) 
p  sinh(pj) 


gm,min)  dp  d^ 


where 


(2.30) 


r 


sm 


kL 


a 

2"^ 


eos  i^x')  dx 


'f 


sm 


kl 


a 

2"^ 


cos(^x)  dx  (2.31) 


and 


COS 


dy' 


cos 

kyfn  2  ) 

dy. 


(2.32) 


In  these  integrals  the  integration  over  spaeial  variables  is  easily  eomputed  using 
trigonometrie  identities  and  integration  of  trigonometrie  funetions.  However,  integrals 
over  speetral  variables  (^,  rj)  are  more  eomplieated.  The  integral  over  rj  can  be  eomputed 
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relatively  easily  using  complex  plane  analysis.  This  process  for  the  example  in  expression 
2.30  is  shown  in  the  next  section. 

2.5  Integrating  ?/  Using  Complex  Plane  Analysis 

In  order  to  evaluate  the  r]  integral  in  expression  2.30,  complex  plane  analysis  is  used. 
Once  gm,mijf)  has  been  evaluated  and  substituted  into  2.30  the  rj  integral  has  the  form 


+(-ir 


cosh(pj)  2  1  -  (-l)'"e 

- 77^ - uT] 

p  sinh(pj) 

cosh(pj)  2  (-l)“-e^''?2Ay 

- 7J  - UTj 

p  sinh(pj)  [jf  _ 


(2.33) 


In  order  to  evaluate  expression  2.33  it  is  first  broken  up  into  5  cases.  These  cases 
involve  differing  values  of  m  and  m.  As  these  values  change,  pole  locations  for  the  function 
will  change,  and  alter  the  outcome  of  the  integral.  The  order  of  the  poles  will  also  change 
as  m  and  m  vary.  Pole  plots  for  these  cases  is  shown  in  Figure  2.5.  Note  that  dp  next  to  a 
pole  indicates  a  double  pole. 


Setting  the  denominator  of  2.33  equal  to  0  identifies  all  singularities.  Poles  77^  and 
p~p  arise  from  the  parallel  plate  structure  due  to  the  Green’s  function  term  sinh(pj)  =  0. 
Note  that  because  of  the  cyclical  nature  of  the  sinh  function,  in  actuality  there  are  an 
infinite  number  of  these  poles.  However,  only  one  is  shown  in  the  pole  plot  for  simplicity. 
Examination  of  the  denominator  of  2.33  also  reveals  there  is  a  double  pole  at  77°  =  0  for 
the  case  of  m  =  m  =  0  (however  it  will  be  shown  that  one  of  these  poles  is  removable). 
In  addition  there  are  either  single  or  double  poles  located  at  77I  =  rjl  =  ^2  -  X’ 

and  J]2  -  depending  on  the  values  of  m  and  in.  There  is  a  branch  point  and  associated 
branch  cut  resulting  from  p  =  0  where  p  is  defined  in  equation  2.27,  but  due  to  the  evenness 
of  the  function  in  p,  it  is  removable  and  therefore  does  not  contribute.  Here,  details  for  the 
evaluation  of  2.33  using  complex  plane  analysis  is  shown  for  the  case  m  =  in  =  0.  Other 
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Figure  2.5:  Pole  locations  for  5  cases  considered  in  the  complex  plane  analysis. 


cases  are  omitted  for  brevity  but  can  be  treated  in  a  similar  manner.  If  m  and  m  are  set  to  0, 
expression  2.33  becomes 


cosh(pd) 
p  sinh(p  J) 

cosh(pj)  I  -  ^ 

p  sinh(p  J) 


(2.34) 


Here  we  examine  the  second  integral  portion  of  2.34.  Cauchy’s  Integral  Theorem 
states  that 
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(2.35) 


f(z)dz  =  0 

for  a  simple  closed  contour  as  long  as  the  function  f(z)  is  analytic  inside  and  on  the  contour. 
Thus  we  are  allowed  to  define  a  contour  around  the  poles  in  figure  2.5.  Closure  in  the  upper 
half  plane  ensures  that  the  conditions  of  Jordan’s  Lemma  [1]  are  satisfied.  This  contour  is 
shown  in  Figure  2.6. 


Figure  2.6:  Contour  for  upper  half  plane  closure  where  m  and  m  equal  0. 


Choosing  this  contour  allows  us  to  use  Cauchy’s  Integral  Theorem  such  that 

d)  f{z)  dz  =  j  +1  +^(p  I  ~  (2.36) 

Jc  J-R  Je  ^  Jc„+  Jc+  Jet 


Using  Cauchy’s  Integral  Formula  for  single  poles,  the  parallel-plate  mode  contribution  from 
the  third  term  on  the  right  hand  side  of  2.36  becomes 


jin  Y  1  - 


(2.37) 
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Since  the  pole  at  the  origin  is  a  double  pole,  Cauchy’s  Integral  Formula  for  higher  order 
poles  is  used  namely, 


fin) 

civ-  Vo)’ 


dr]  = 


jn  d"  V(7/) 


{n  -  1)!  di] 


n-\ 


ri=m 


for  a  half  circle  where  n  is  the  order  of  the  pole.  In  this  problem 


(2.38) 


fiv) 

and  n  =  2.  Using  2.38  on  2.39  yields 
cosh(pj) 


coshipd) 
p  sinhipd) 


(2.39) 


p  sinh(pd) 


d  cosh(pd) 


77=0 


di]  p  sinh(jc>  J) 


77=0  J 


(2.40) 


It  is  easy  to  see  that  the  second  half  of  expression  2.40  goes  to  0  when  it  is  evaluated  at 


p  =  Po  =  0.  Thus  the  contribution  from  Ct  becomes 


2nAy 


coshipod) 


(2.41) 


Po  sinhipod) 

where  po  =  -  k^.  It  is  also  easily  shown  using  the  substitution  p  =  that  the 

integrand  from  adheres  to  Jordan’s  Lemma  and  therefore  has  0  contribution  as  ^  oo. 
Finally  using  Cauchy’s  Integral  Theorem  in  the  limit  as  R  approaches  oo  and  5  approaches 
0"^,  it  can  be  shown  that 


lim  lim 

^—>0+  R-^oo 


X— (5  ^oo 

+  j  =  PV  =  27TAy- 
R  JS  j  J-oo  i 


coshipod)  j2n 


--tL  '  (2.42) 


Po  smh(pod) 

where  PV  denotes  the  integral  has  been  computed  in  a  principle  value  sense.  Note  that  as 
R  approaches  oo  more  parallel  plate  poles  are  included.  The  first  half  of  2.34  uses  the  same 
process  as  the  second  half.  However,  note  that  this  portion  must  be  closed  in  the  lower  half 
plane  to  ensure  satisfaction  of  Jordan’s  Lemma.  A  contour  for  lower  half  closure  is  shown 
in  Figure  2.7 

After  using  the  process  detailed  above  it  can  be  shown  that  the  first  half  of  2.34 
evaluates  to  the  same  expression  as  the  second  half  of  2.34.  Thus  the  final  form  of  2.34  is 
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Figure  2.7:  Contour  for  lower  half  plane  closure  where  m  and  m  equal  0. 


simply  twice  the  value  of  equation  2.42.  Once  this  case  is  complete,  the  4  other  cases  for 
this  example  are  evaluated  using  the  same  steps.  The  63  other  integrals  are  also  evaluated 
using  the  same  steps. 

2.6  Calculating  F  and  T 

Once  all  t]  integrals  are  evaluated,  the  ^  integrals  are  computed  numerically.  These 
integrals  contain  branch  cuts  and  cannot  be  found  in  closed  form.  MATLAB  was  used  for 
all  code  and  simulations  for  this  research.  Gaussian  quadrature  was  used  to  integrate  all  ^ 
integrals.  Once  all  ^  integrals  are  calculated  a  linear  system  of  equations  remains  where  F 
and  T  are  the  unknowns.  This  system  can  be  rearranged  such  that: 


<3l,l 

(^2,\ 


<^\,2 

Cl2,2 


bl 

b2 


(2.43) 
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2.7  Simulated  Dual  Probe  DRWG  Results 


Results  for  simulating  the  dual  probe  DRWG  system  are  shown  below.  Fixed 
parameters  of  the  system  are  shown  in  Table  (2.1).  The  MUT  for  this  simulation  was 
FGM125,  a  thin,  flexible,  broadband,  magnetieally  loaded  silieone  absorber  [28]. 


Table  2.1:  Simulation  parameters  for  the  dual  probe  DRWG  system  showing  dimensions 
for  the  DRWGs  and  fgml25  parameters. 


a 

.0183  m 

b 

.0082  m 

Av: 

.0022  m 

Ay 

.0013  m 

Jo 

(1.01  xb)m 

d 

.0031  m 

MUT 

fgml25 

U 

6.7314-0.22251 

fir 

0.6202  -  0.2960i 

f 

6.0  GHz 

Initial  simulation  results  were  used  to  eonfirm  the  MATLAB  eode  is  working.  A  result 
for  F  and  T  was  obtained  and  was  physieally  realizable.  It  was  expeeted  that  the  magnitude 
of  both  refleetion  and  transmission  eoeflieients  would  range  somewhere  between  0  and  1 . 
After  the  initial  test  eoneluded,  another  series  of  simulations  were  eondueted  to  further  test 
the  system.  The  first  of  these  tests  eonsisted  of  varying  yo-  Beeause  of  the  lossy  material 
simulated,  the  expeetation  is  that  by  inereasing  y^,  T  should  deerease;  and  by  deereasing 
yo,  T  should  inerease.  These  results  are  shown  in  table  (2.2).  The  seeond  test  eonsisted  of 
varying  fir-  The  expeetation  is  that  by  deereasing  the  imaginary  eomponent  of  fir  (making 
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the  material  low  loss)  T  will  increase  because  less  energy  is  lost  in  the  MUT.  These  results 
are  shown  in  table  (2.3). 


Table  2.2:  Results  of  Varying  yo  showing  decreases  in  T  as  yo  increases  for  lossy  media 


yo 

\T\ 

1.0001  xZ? 

.0502 

imxb 

.0499 

1.5  xb 

.0367 

\0xb 

1.9 

Table  2.3:  Results  of  Varying  showing  increases  in  T  as  fir  decreases  for  lossy  media 


fir 

in 

0.6202  -  0.2960/ 

.0499 

0.6202-0.1/ 

.0659 

0.6202  -  0.01/ 

.0779 

Tables  2.2  and  2.3  show  all  results  are  trending  as  expected.  It  also  shows  that,  while 
using  lossy  media,  T  decreases  to  essentially  0  if  the  probes  are  placed  too  far  apart.  This 
implies  that  when  using  lossy  media  the  two  DRWG  probes  need  to  be  as  close  as  possible 
to  ensure  sufficient  coupling  to  obtain  an  accurate  T  measurement.  A  further  simulation 
was  conducted  to  determine  if  the  minimum  physical  distance  between  waveguides  is 
sufficient  to  measure  T  in  lossy  media.  In  this  simulation  the  same  lossy  MUT  was  used 
at  a  frequency  of  12  GHz.  Waveguide  offset  yo  was  then  varied  in  finer  increments.  These 
results  are  shown  in  Figure  (2.8). 
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Figure  2.8:  F  and  T  versus  3;o  showing  the  minimum  distanee  between  waveguides  for  lossy 
material 


These  simulation  results  eonfirm  that  the  minimum  realizable  yo  distanee  denoted  by 
the  red  lines  will  ensure  a  measurable  transmission  eoeffieient  for  lossy  media.  It  also 
shows  oseillating  behavior  in  the  refleetion  measurement.  This  ean  be  attributed  to  eoupling 
between  the  two  waveguides  going  in  and  out  of  phase.  A  final  simulation  was  eondueted 
to  determine  how  the  system  would  behave  when  a  low  loss  material  was  used  instead  of 
a  high  loss  material.  The  low  loss  material  simulated  was  Teflon  with  material  parameters 
6r  =  2.1  -  y.OOl  and  =  T  The  simulation  used  the  same  waveguide  parameters  at  a 
frequeney  of  12  GHz.  Results  of  this  simulation  are  shown  in  Figure  2.9. 

Figure  2.9  shows  that,  as  expeeted,  T  does  not  deerease  as  fast  for  low  loss  material. 
This  allows  for  greater  spaeing  between  waveguides  when  taking  low  loss  measurements.  It 
may  also  be  benefieial  to  have  the  probes  farther  apart  to  allow  the  waves  to  interaet  further 
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Figure  2.9:  F  and  T  versus  jo  for  low  loss  material. 


in  distanee  with  the  low  loss  material.  Again  oseillation  is  seen  in  the  measurements,  a 
further  sign  of  eonstruetive  and  destruetive  interferenee  eoupling  between  the  waveguides. 
Beeause  of  this  oseillation,  it  makes  determining  the  best  spaeing  between  probes  more 
difficult  for  low  loss  material.  It  is  recommended  that  the  probes  be  spaced  in  such  a 
way  that  the  measurement  occurs  on  the  flatter  part  of  the  curves.  This  will  help  reduce 
measurement  errors  associated  with  errors  in  yo-  However,  the  flatter  part  of  the  curve 
will  move  as  the  frequency  changes.  This  may  render  choosing  the  spacing  based  on 
minimization  of  errors  due  to  yo  more  difficult. 

Comparing  Figures  2.8  and  2.9,  both  show  exponential  decay  in  the  transmission 
measurement.  Differences  in  decay  rate  are  driven  by  the  amount  of  loss  and  wave 
spreading  in  each  material.  However,  in  the  event  of  a  lossless  material,  there  would  still 
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be  algebraic  decay  based  on  wave  spreading.  Upon  consideration  of  the  findings  in  Figures 
2.8  and  2.9,  it  is  concluded  that  sufficient  coupling  does  occur  for  both  lossy  and  lowloss 
materials.  Thus,  the  dual  probe  DRWG  comprises  a  valuable  candidate  as  a  fiight-line 
diagnostic  tool. 
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III.  Multimode  Single  Probe  DRWG 


In  this  chapter,  the  analysis  of  a  single  probe  DRWG,  shown  in  Figure  (3.1),  is 
performed.  The  analysis  includes  an  examination  of  the  reflection  coefficients  of  the 
dual  probe  DRWG  placed  on  a  conductor-backed  material.  It  is  assumed  in  the  analysis 
that  only  a  single  dominant  mode  is  propagating  in  the  free-space  filled  DRWG  however, 
higher  order  evanescent  modes  are  allowed  to  exist  during  the  analysis.  The  fundamental 
goal  for  performing  this  analysis  is  to  assess  if  a  single  mode  is  sufficient  to  accurately 
determine  permittivity  and  permeability  of  the  material,  or  whether  higher  order  modes 
must  be  included  in  the  analysis  to  obtain  accurate  permittivity  and  permeability  results. 
This  answer  will  greatly  impact  the  ability  to  perform  these  measurements  in  real  time, 
versus  needing  a  large  amount  of  time  and  resources  to  post-process. 

The  multimode  single  probe  DRWG  analysis  uses  the  same  field  expressions  detailed 
in  chapter  2.  First  field  expressions  in  the  DRWG  are  defined.  Next,  the  fields  in  the  PP 
region  of  the  single  probe  DRWG  are  found  with  the  aid  of  Love’s  Equivalence  Theorem. 
Subsequent  enforcement  of  tangential  boundary  conditions  at  the  single  probe  DRWG 
aperture  leads  to  a  MFIE  for  the  unknown  reflection  coefficient.  Recall  however,  that  if 
two  material  properties  are  desired,  two  linearly-independent  measurements  are  required. 
To  obtain  this  second  independent  measurement,  a  two  thickness  technique  is  used  where 
another  reflection  measurement  is  taken  for  the  same  material  of  a  different  thickness. 
The  two  coupled  MEIE’s  are  solved  via  the  Method-of-Moments  for  a  lossy  material. 
Experimental  results  for  the  calculation  of  6^  and  /Ur  are  presented  and  discussed,  and 
general  conclusions  are  provided. 
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Figure  3.1:  Side  view  of  single  probe  DRWG  geometry 


3.1  Single  Probe  DRWG  Region  Fields 

As  a  wave  travels  down  the  DRWG  and  internets  with  the  boundary  between  the 
waveguide  and  the  MUT,  a  portion  of  the  energy  will  be  refleeted  baek  into  the  waveguide 
[5].  Thus  we  must  first  find  the  total  fields  in  the  DRWG  region.  All  fields  in  a  single  probe 
DRWG  system  are  found  using  a  proeess  similar  to  that  used  in  ehapter  2.  Total  fields  in 
the  DRWG  portion  of  the  system  are  found  using  the  DRWG  fields  developed  in  ehapter  2. 
While  ehapter  2  only  uses  a  single  propagating  mode,  a  more  robust  multi-mode  analysis 
allows  for  multiple  refleetion  eoeffieients  for  many  potential  waveguide  modes.  Allowing 
for  this  yields  fields  in  the  form  of 


where 


N 

'Edrwg  =  ri(p )  +  ^  ,q(p) 

?=1 

N 

'Hdrwg  =  e~^''-^^~haCp)  -  ^  Yqe^’^^^'h  tq(p) 

?=i 


e  t  =  xE,,  +  yEy 


(3.1) 

(3.2) 


(3.3) 
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(3.4) 


Variable  is  the  reflection  coefficient  associated  with  the  even  mode.  Only  even 
modes  are  considered  because  the  dominant  mode  is  even  symmetric  and  the  system  is  even 
symmetric.  This  forces  all  possible  excited  modes  to  be  even.  Variable  N  is  the  number 
of  modes  allowed  to  be  reflected  and  can  be  as  small  or  large  as  required  for  convergence. 
Upon  inspection  it  is  easy  to  see  that  equations  3.1  and  3.2  simplify  to  an  equation  similar 
to  that  of  2.6  and  2.7  if  V  =  1  and  only  one  reflected  mode  is  considered. 


3.2  Parallel  Plate  Region  Fields 

As  a  wave  travels  down  the  DRWG  and  contacts  the  interface  between  the  waveguide 
and  the  MUT,  a  portion  of  the  energy  will  be  transmitted  into  the  MUT.  This  energy  gives 
rise  to  electromagnetic  fields  inside  the  PP  region  of  the  system.  These  fields  are  found 
using  a  process  similar  to  that  used  in  chapter  2.  First  Love’s  equivalence  principle  [5]  is 
used  to  replace  the  actual  source  (the  DRWG  aperture)  with  an  equivalent  source  in  the 
form  of  magnetic  surface  currents,  namely  [5] 


M  =  -hxE.  (3.5) 

where  h  the  normal  unit  vector  points  into  the  MUT  region.  This  equivalent  system  is 
shown  in  Figure  (4.2). 


Now  that  the  equivalent  current  has  been  identified,  a  vector  potential  (VP)  solution  is 
used  to  find  the  parallel  plate  region  fields.  This  solution  is  as  follows: 

t  +  V,V-  )FP/^,z).  (3.6) 

joje/j. 

However,  the  F  VP  must  first  be  identified.  The  surface  current  Mi  supports  this  vector 
through  the  following  relation. 


Ff^(p,z)=  rr  G(p,z[p' ,0)- eMi(p')ds' 


kJ  U  s  I 


(3.7) 
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Figure  3.2:  Equivalent  system  for  parallel  plate  region  of  the  single  probe  DRWG  system 
after  Love’s  Equivalenee  Theorem  is  applied. 


where  G0 ,z[p  ' ,0)  is  the  parallel  plate  Green’s  funetion,  and  Mi  is  found  using  equation 
3.5.  Equating  the  tangential  magnetie  field  at  the  aperture  leads  to  the  MEIE 

- - (k^I  t  +  VtV-jFf‘’(p,z)  =  ha(p)-'y^T'ghtg(p).  (3.8) 

This  MEIE  is  solved  in  the  next  seetion  via  the  MoM. 

3.3  Solving  the  MFIE 

This  seetion  details  solving  the  MEIE  of  the  previous  seetion.  The  first  step  of  the 

— ^ 

MoM  is  to  expand  the  unknown  eurrent  (Mj)  as  a  series  of  basis  funetions.  In  this  ease 
we  will  exploit  the  faet  that  we  know  the  fields  in  the  DRWG  region  just  to  the  left  of  the 
boundary,  and  we  know  that  these  fields  must  be  eontinuous  aeross  the  boundary,  namely 

MiCp  )  =  -h  X  E  X  -h  X  Edrwg-  (3.9) 

Substituting  3.9  into  3.7  and  further  substituting  into  3.8  yields  the  following  relation 
after  expansion: 
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+ V(V-)  JJ'  G(p,  zfp', 0)- X  EoRwcjds'  =  htiG))-^r^htq(p)  (3.10) 

Next  both  sides  of  equation  3.10  are  tested  with  N  total  magnetic  field  modes.  This 
must  be  done  N  times  because  there  are  N  unknowns  in  the  form  of  Ti  through  T^?.  The 
testing  operator  is  as  follows: 

ffs)- 

After  testing,  3.10  becomes 


h*dS, 


(3.11) 


JJ  +  JJ  G(p,z[p',0)-(^-nx~EoRWGjds'j-~h*gdS  I 

=  11  ha  - 


(3.12) 


h*dS,. 


From  this  point  forward  it  is  desirable  to  work  with  the  spectral  representation  of  the 
Green’s  function.  The  Green’s  function  is  shown  below 


G0,z[p',O)  =  ^ff_  G('A,zl0)e^^^^-^'0eJ^<^-y’0d^drj 


where 


and 


G(d,z|0)  =  ^ 


cosh(pj) 


p  sinh(pd) 

p  =  +  rf  -  k^. 


(3.13) 


(3.14) 


(3.15) 


After  using  this  Green’s  function,  rearranging  the  order  of  integration,  and  pulling  out 
constants,  equation  3.12  becomes 


‘  rr  rr  rr 

jcop(2ny  JJ_oo  JJs,  JJs, 


(k  7  I  +  VfV  ■)  G(A,  z|0)  •  -hx  E  drwg  ) )  •  h*^ds'  ds  drj  d^ 


(3.16) 


Once  all  operations  are  conducted  on  3.16,  a  series  of  six  dimensional  integrals  remain 
on  the  left  of  the  equality,  and  a  double  integral  remains  on  the  right  of  the  equality.  The 
double  integral  on  the  right  has  three  sections  that  must  be  calculated  separately,  one  for 
each  region  of  the  waveguide.  These  integrals  are  straightforward  and  easily  calculated 
using  standard  trigonometric  identities  and  integration  of  trigonometric  functions.  Note 
these  integrals  must  be  calculated  N  times  due  to  there  being  N  unknowns.  Because  the 
waveguide  modes  are  orthogonal  only  values  where  expansion  and  testing  indices  are  equal 
(q  =  q)  need  to  be  calculated.  The  following  expression  is  the  final  solution  to  the  right 
side  double  integral  of  equation  3.16  where  ^  =  ^  =  1. 


l-Ti 


IZrsp 


Xi  (sinc(2jlf„  ,.^A.r)  +  sinc(2fcf„  ,,A.r))  AxAy(l  +  d„,o) 

n 

Yj  {smc(2jkl^  .^Ax)  -  sinc(2fcf„  „A.r))  AxAy(l  -  d„,o) 

n 

+  Xi  Zi  -  A.r 


n  n,m  ' 

-sine  ( 2fc^^ -  A;c))  j  -  A.r)  ^(1  +  6^^) 


+  Z  Z  sine  (t  “ 


(3.17) 


+sinc  ( 2fc(^ ,,  -  A.r))  ) 


The  left  side  integrals  are  more  complicated.  In  the  single  probe  system  16  total  6 
dimensional  integrals  must  be  calculated  for  each  combination  of  q  and  q  leading  to  a  total 
of  77^  X 16  six  dimensional  integrals.  One  example  of  these  6  dimensional  integrals  is  shown 
below 


ZZ 

n,m  h,m 


t  t* 

^yn,m  ^yn,m 


Xoo 

OO  %J —0( 


cosh(pj) 
p  sinh(pj) 


gm,miv)  dq 


(3.18) 
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where 


f 


Sin 


kL 


--X' 
2  "" 


cos(^y)  dx‘ 


Sin 


kl, 


a 

2-\ 


cos,{^x)  dx  (3.19) 


and 


(  ,  b\ 

cos 

kym  “  2 ) 

dy’ 


cos 


^ym 


dy. 


(3.20) 


In  this  example  integrals  with  respect  to  spatial  variables  (x,x',y,y')  are  easily 
computed  using  trigonometric  identities  and  integration  of  trigonometric  functions. 
However,  integrals  over  spectral  variables  (^,  r])  are  more  complicated.  The  integral  over 
rj  can  be  computed  relatively  easily  using  complex  plane  analysis  i.e.  Cauchy’s  theorem. 
For  more  information  on  complex  plane  integration  see  chapter  2.  Once  all  rj  integrals  are 
computed  the  ^  integrals  are  computed  numerically.  These  integrals  contain  branch  cuts 
and  cannot  be  found  in  closed  form.  MATLAB  was  used  for  all  code  and  simulations  for 
this  research.  Gaussian  quadrature  was  used  to  integrate  all  ^  integrals.  Once  all  ^  integrals 
are  calculated  a  linear  system  of  equations  remains  where  T i  through  T n  are  the  unknowns. 
This  system  can  be  rearranged  such  that: 
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(3.21) 


It  is  important  to  note  that  although  there  are  multiple  unknowns  in  equation  3.21, 
these  values  are  not  typically  measurable.  Only  T i  is  typically  measured  since  it  represents 
reflection  of  the  dominant  propagating  mode.  In  order  to  get  a  second  independent 
measurement  and  solve  for  both  permittivity  and  permeability,  a  two  thickness  technique  is 
used  to  get  the  second  measurement.  This  second  measurement  uses  the  exact  same  process 
as  previously  described  where  the  only  difference  is  the  thickness  of  the  MUT.  Because  the 


36 


derivation  is  exactly  the  same,  no  new  theoretical  calculations  need  to  be  done.  The  exact 
same  MATLAB  code  used  to  determine  equation  3.21  is  used  with  a  second  thickness. 
This  second  thickness  for  this  research  is  simply  double  the  original  thickness.  To  take 
experimental  measurements,  two  layers  of  MUT  are  used  instead  of  one.  This  second 
linear  system  of  equations  is  arranged  such  that: 
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where  superscript  2  denote  computations  for  thickness  2.  Once  equations  3.21  and  3.22  is 
complete  material  parameters  e  and  can  be  calculated. 


3.4  Calculating  e  and  fx 

Now  that  all  calculations  needed  to  obtain  and  solve  equations  3.21  and  3.22  are 
complete,  material  parameters  e  and  jj.  can  be  computed.  In  this  step  these  parameters 
are  solved  iteratively  using  a  least  squares  technique.  First  an  initial  guess  for  e  and  //  is 
used  to  calculate  initial  T}  and  T^  estimates.  These  estimates  are  compared  to  experimental 
measurements  taken  from  a  single  probe  system.  Experimental  measurement  data  used 
was  the  same  data  obtained  and  used  in  [21].  After  comparing  T}  and  Tj  to  reflection 
coefficients  of  the  experimental  system,  an  update  to  e  and  jj.  is  calculated  and  used  to 
recalculate  Tj  and  Tj.  This  process  happens  repeatedly  until  the  difference  in  Tj  and  T^ 
and  the  experimental  data  is  within  tolerances  set  by  the  user.  MATLAB ’s  Isqcurvefit 
function  was  used  to  perform  the  iterative  least  squares  operation.  Experimental  results  of 
this  process  are  shown  in  the  next  section. 
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3.5  ECCOSORB  SF-3  Experimental  Results 

This  section  details  the  experimental  results  using  multiple  modes  for  the  single  probe 
DRWG  system.  All  experimental  data  came  from  [21],  see  this  for  more  details  on  the 
experimental  setup  and  equipment  used.  In  this  experiment,  material  characterization 
measurements  were  made  at  6-18  GHz  of  a  lossy  silicon-based  magnetic  absorbing  material 
ECCOSORB  SF-3  [16]  manufactured  by  Emerson  &  Cuming  Microwave  Products,  Inc. 
Figures  (3.3)  and  (3.4)  show  calculated  e  and  n  values  respectively.  Both  figures  show 
calculations  for  1  through  3  modes,  and  also  experimental  results  using  a  Nicholson-Ross- 
Weir  (NRW)  method  [22]. 


s  and  s  for  SF-3  Single  Probe 


Figure  3.3:  Calculated  e  values  for  1-3  modes  as  well  as  experimental  results  using  a 
NRW  method. 
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|i^  and  for  SF-3  Single  probe 


Figure  3.4:  Calculated  yu  values  for  1  -  3  modes  as  well  as  experimental  results  using  a 
NRW  method. 


Figures  3.3  and  3.4  show  that  e  and  yU  calculations  using  multiple  modes  are  consistent 
with  those  using  a  single  mode.  They  are  also  consistent  with  the  results  using  the 
NRW  method.  The  small  differences  in  results  can  be  attributed  to  measurement  errors 
and  sample  inconsistency.  While  using  higher  modes  appears  to  yield  a  better  result 
(converges  closer  to  the  NRW  method)  for  some  frequencies  this  is  not  the  case  in  all  of  the 
frequencies.  This  leads  to  the  conclusion  that  adding  higher  order  modes  when  calculating 
6  and  yu  does  not  yield  a  significant  benefit  over  using  a  single  mode.  This  also  proves 
that  for  this  particular  system,  higher  order  modes  are  not  being  excited.  If  they  are  being 
excited,  the  amount  of  energy  coupling  is  extremely  low.  This  conclusion  is  confirmed  by 
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examining  the  expected  reflection  coefflcients  for  each  mode.  These  expected  coefficients 
(calculated  with  equations  3.21  and  3.22)  are  shown  in  Table  3.1,  and  show  that  the  amount 
of  energy  reflected  in  the  second  and  third  mode  are  1  and  2  orders  of  magnitude  lower 
than  that  of  the  dominant  mode. 


Table  3.1:  Expected  reflection  coefficients  for  each  mode 


Mode 

ini 

|r2l 

1 

.3536 

.5265 

2 

.0299 

.0369 

3 

.0026 

.0017 

Upon  consideration  of  the  findings  in  Figures  3.3  and  3.4,  as  well  as  Table  3.1,  it 
is  concluded  that  higher  order  modes  do  not  need  to  be  used  for  this  class  of  materials  to 
calculate  accurate  material  parameters.  Only  needing  a  single  mode  will  save  valuable  time 
and  resources,  and  will  yield  a  valuable  diagnostic  tool. 
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IV.  Multimode  Clamped  Probe  DRWG 


In  this  chapter,  the  analysis  of  a  clamped  probe  DRWG,  shown  in  Figure  4.1,  is 
performed.  The  analysis  includes  an  examination  of  the  reflection  and  transmission 
coefficients  of  a  lossy  material  sandwiched  between  two  DRWGs.  It  is  assumed  in  the 
analysis  that  only  a  single  dominant  mode  is  propagating  in  the  free-space  filled  DRWG’s 
however,  higher  order  evanescent  modes  are  allowed  to  exist  during  the  analysis.  The 
fundamental  goal  for  performing  this  analysis  is  to  assess  if  a  single  mode  is  sufficient  to 
accurately  determine  permittivity  and  permeability  of  the  material,  or  whether  higher  order 
modes  must  be  included  in  the  analysis  to  obtain  accurate  permittivity  and  permeability 
results.  This  answer  will  greatly  impact  the  ability  to  perform  these  measurements  in  real 
time,  versus  needing  a  large  amount  of  time  to  post-process. 

The  multimode  clamped  probe  DRWG  analysis  uses  the  same  field  expressions 
detailed  in  chapter  2.  First  field  expressions  in  the  DRWG’s  are  defined.  Next,  the  fields  in 
the  PP  region  of  the  clamped  probe  DRWG  are  found  with  the  aid  of  Love’s  Equivalence 
Theorem.  Subsequent  enforcement  of  tangential  boundary  conditions  at  the  clamped  probe 
DRWG  apertures  leads  to  coupled  MFIE’s  for  the  unknown  reflection  and  transmission 
coefficients.  These  MEIE’s  are  solved  via  the  Method-of-Moments  for  a  lossy  material. 
Experimental  results  for  the  calculation  of  6^  and  fir  are  presented  and  discussed,  and 
general  conclusions  are  provided. 

4.1  Clamped  Probe  DRWG  Region  Fields 

As  a  wave  travels  down  the  first  DRWG  and  interacts  with  the  boundary  between  the 
waveguide  and  the  MUT,  a  portion  of  the  energy  will  be  reflected  back  into  the  waveguide 
[5].  Energy  will  also  be  transmitted  through  the  PP  region  and  into  the  second  waveguide. 
Thus  we  must  first  find  the  total  fields  in  the  DRWG  regions.  All  fields  in  a  clamped  probe 
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Figure  4.1:  Side  view  of  clamped  probe  DRWG  geometry 


DRWG  system  are  found  using  a  process  similar  to  that  used  in  chapter  2.  Total  fields  in 
the  DRWG  portions  of  the  system  are  found  using  the  DRWG  fields  developed  in  chapter 
2.  While  chapter  2  only  uses  a  single  mode,  a  more  robust  multi-mode  analysis  allows 
for  multiple  reflection  and  transmission  coefficients  for  many  potential  waveguide  modes. 
Allowing  for  this  yields  fields  in  the  form  of 


where 


N 

Edrwg\  =  ti(p)  +  ^  tq(p) 

q=\ 

N 

'HdRWGI  =  ^  fl(p  )  -  ^  h  tq(p) 

q=l 

N 

EdRWGI  =  tqCp) 

q=l 
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Hdrwgi  =  J]Tqe~^'^-^^^-‘^'>'hrq(p) 

q=l 


e  t  =  xE^  +  yEy 


zx  e  t 
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(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 
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Variable  is  the  reflection  coefficient  associated  with  the  even  mode,  and  is  the 
transmission  coefficient  associated  with  the  even  mode.  Only  even  modes  are  considered 
because  the  dominant  mode  is  even  symmetric  and  the  system  is  even  symmetric.  This 
forces  all  possible  excited  modes  to  be  even.  Variable  N  is  the  number  of  modes  allowed  to 
be  reflected  and  transmitted  and  can  be  as  small  or  large  as  required  for  convergence.  Upon 
inspection  it  is  easy  to  see  that  equations  4.1  and  4.2  simplify  to  an  equation  similar  to  that 
of  2.6  and  2.7  if  V  =  1  and  only  one  reflected  mode  is  considered. 

4.2  Parallel  Plate  Region  Fields 

As  a  wave  travels  down  the  first  DRWG  and  contacts  the  interface  between  the 
waveguide  and  the  MUT,  a  portion  of  the  energy  will  be  transmitted  into  the  MUT.  This 
energy  gives  rise  to  electromagnetic  fields  inside  the  PP  region  of  the  system.  These 
fields  are  found  using  a  process  similar  to  that  used  in  chapter  2.  First  Love’s  equivalence 
principle  [5]  is  used  to  replace  the  actual  sources  (the  DRWG  apertures)  with  equivalent 
sources  in  the  form  of  magnetic  surface  currents,  namely  [5] 

Mi=-flx'Ei.  (4.7) 

M2  = -fix  £2.  (4.8) 

where  h  the  normal  unit  vector  points  into  the  MUT  region.  Note  that  h  points  in  the 
positive  z  direction  for  4.7  and  the  negative  z  direction  for  4.8.  This  equivalent  system  is 
shown  in  Figure  (4.2). 

Now  that  the  equivalent  currents  have  been  identified,  a  VP  solution  is  used  to  find  the 
parallel  plate  region  fields.  This  solution  is  as  follows: 

MrG^,z)  =  —(k^I  t  +  V,V-  )FrG^,z).  (4.9) 

jcoeix 
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Figure  4.2:  Equivalent  system  for  parallel  plate  region  of  the  elamped  probe  DRWG  system 
after  Love’s  Equivalenee  Theorem  is  applied. 


However,  the  F  VP  must  first  be  identified.  The  surfaee  eurrents  M\  and  M2  support  this 
veetor  through  the  following  relation, 


Ff^Go,z)=  ff  G^(p,z[p',0)-eMi(p')ds'  +  ff  ,z[p' ,d)- 6M2(^')ds'  (4.10) 

JJSl  JJs2 


where  G^Cp,  z\p  ' ,G)  and  G^(^ ,  z[p' ,d)  are  the  parallel  plate  Green’s  funetions  assoeiated 
with  Ml  and  M2  respeetively,  and  Mi  and  M2  are  found  using  equation  4.7  and  4.8. 
Equating  the  tangential  magnetie  field  at  the  apertures  leads  to  the  MEIEs 


jojefi 


N 

VfV-)Ff^(p,z)  =  hn(p)-^T'ghtq(p). 

g=l 


—  (k^It  + 
jcoen 


N 

g=l 


These  MEIEs  are  solved  in  the  next  seetion  via  the  MoM. 


(4.11) 

(4.12) 


4.3  Solving  the  MFIEs 

This  seetion  details  solving  the  MEIEs  of  the  previous  seetion.  The  first  step  of  the 
MoM  is  to  expand  the  unknown  eurrents  (Mi)  and  (M2)  as  a  series  of  basis  funetions.  In 
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this  case  we  will  exploit  the  fact  that  we  know  the  fields  in  the  DRWG  region  just  to  the 
outside  of  the  boundaries,  and  we  know  that  these  fields  must  be  continuous  across  the 
boundaries,  namely 


)  —  —h  X  E I  ^  —h  X  Ej)]iWG\-  (4.13) 

M2(p')  =  -hx^2~  -ux'Edrwgi-  (4.14) 

Substituting  4.13  and  4.14  into  4.10  and  further  substituting  into  4.11  and  4.12  yields  the 
following  relations  after  expansion: 


(k^/,  +  V,V-)  rr  G^(p,z\p' ,0)-Mi(p')ds'  +  ff  G2(p,z[p',d)-M2(p')ds' 
JJsi  JJs2 


N 


hn(p) 

Cl=l 

(4. 


-(k2/,  +  V,V-) 


rr  G^(p,z[p' ,0)-Mi(p')ds'  +  rr  G^2(p,z[p' ,d)-M2(p')ds' 
jjsi  JJs2 


5) 


N 


=  J^TghtqCp) 


g=l 


(4.16) 


Next  both  sides  of  equations  4.15  and  4.16  are  tested  with  N  total  magnetic  field 
modes.  This  must  be  done  N  times  because  there  are  N  unknowns  in  each  equation  in  the 
form  of  Ti  through  T^  and  Ti  through  The  testing  operators  are  as  follows: 


J£{---}-7^:qdS2 

for  4.15  and  4.16  respectively.  After  testing  4.15  and  4.16  become 


(4.17) 

(4.18) 
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rr  + v,v-)  rr  G^(^,z\p\())-Mi(^')ds' 

JjsAJ^d  [JJsi 

+  JJ'  G2(p,z[p',d)-M2(p')ds'  j-  h*^dSi 

^  jj  ~  X  1 

JJ^  G^(p,z[p',0)-Mi(p')ds' 


T  I  ^ 

5  2  \ 


+ 


h%dS2 


JJ'  G2(p,zfp',d)-M2(p')ds' 

=  JJ  (ET,X,(7^)j-t;,J5 


(4.19) 


(4.20) 


From  this  point  forward  it  is  desirable  to  work  with  a  spectral  representation  of  the 
Green’s  function.  The  Green’s  functions  are  shown  below 


where 


and 


G^(^,z[p',0)  =  ^ff_ 
G,(^,z[p',d)  = 


Gi(d,z|0)  =  /, 

G,(A,zld)  =  [, 


cosh(pd) 
p  sinh(p  J) 
1 

p  sinh(p  J) 


p  =  +  rf  -  k^. 


(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 
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After  using  this  Green’s  functions,  rearranging  the  order  of  integration,  and  pulling 


out  constants,  equations  4.19  and  4.20  become 

‘  rr  ff 

jojuiln)^  JJ_oo  JJs, 
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(4.27) 


Ills  Tjh,g(p)j-7i*^dS: 


Once  all  operations  are  conducted  on  4.26  and  4.27,  a  series  of  six  dimensional 
integrals  remain  on  the  left  of  the  equalities,  and  a  double  integral  remains  on  the  right 
of  the  equalities.  The  double  integral  on  the  right  has  three  sections  that  must  be  calculated 
separately,  one  for  each  region  of  the  waveguide.  These  integrals  are  straightforward  and 
easily  calculated  using  standard  trigonometric  identities  and  integration  of  trigonometric 
functions.  Note  these  integrals  must  be  calculated  N  times  due  to  there  being  N  unknowns. 
Because  the  waveguide  modes  are  orthogonal  only  values  where  expansion  and  testing 
indices  are  equal  (q  =  q)  need  to  be  calculated.  The  left  side  integrals  are  more  complicated. 
In  the  clamped  probe  system  64  total  6  dimensional  integrals  must  be  calculated  for  each 


47 


combination  of  q  and  q  leading  to  a  total  of  64  six  dimensional  integrals.  An  example 
of  one  of  these  integrals  can  be  seen  in  chapters  2  and  3. 

In  these  integrals  the  integration  over  spacial  variables  {x,  x',y,  y')  are  easily  computed 
using  trigonometric  identities  and  integration  of  trigonometric  functions.  However, 
integrals  over  spectral  variables  (^,  q)  are  more  complicated.  The  integral  over  q  can  be 
computed  relatively  easily  using  complex  plane  analysis  (i.e.  Cauchy’s  theorem).  For 
more  information  on  complex  plane  integration  see  chapter  2.  Once  all  q  integrals  are 
computed  the  ^  integrals  are  computed  numerically.  These  integrals  contain  branch  cuts 
and  cannot  be  found  in  closed  form.  MATLAB  was  used  for  all  code  and  simulations  for 
this  research.  Gaussian  quadrature  was  used  to  integrate  all  ^  integrals.  Once  all  ^  integrals 
are  calculated  a  linear  system  of  equations  remains  where  Ti  through  T^r  and  Tj  through 
Ta?  are  the  unknowns.  This  system  can  be  rearranged  such  that: 


«i,i 


Ti 

bi 

Tn 

Ti 

T;v 

b2N 

(4.28) 


Y<^2N,\  .  <^2N,2n\ 

It  is  important  to  note  that  although  there  are  multiple  unknowns  in  equation  4.28, 
these  values  are  not  typically  measurable.  Only  Ti  and  Tj  are  typically  measured  since 
they  represent  reflection  and  transmission  of  the  dominant  propagating  mode.  Unlike  the 
single  probe,  no  more  calculations  need  to  be  made  as  Ti  and  Ti  are  two  independent 
measurements. 
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4.4  Calculating  e  and  ju 

Now  that  all  calculations  needed  to  obtain  and  solve  equation  4.28  are  complete, 
material  parameters  e  and  jj.  can  be  computed.  In  this  step  these  parameters  are  solved 
iteratively  using  a  least  squares  technique.  First  an  initial  guess  for  e  and  //  is  used 
to  calculate  initial  Fi  and  Ti  estimates.  These  estimates  are  compared  to  experimental 
measurements  taken  from  a  clamped  probe  system.  Experimental  measurement  data  used 
was  the  same  data  obtained  and  used  in  [20].  After  comparing  Fi  and  Ti  to  reflection  and 
transmission  coefficients  of  the  experimental  system,  an  update  to  e  and  /i  is  calculated  and 
used  to  recalculate  F i  and  Ti .  This  process  happens  repeatedly  until  the  difference  in  Fj  and 
Ti  and  the  experimental  data  is  within  tolerances  set  by  the  user.  Broy den’s  Method  was 
used  to  perform  the  iterative  least  squares  operation.  Experimental  results  of  this  process 
are  shown  in  the  next  section. 

4.5  ECCOSORB  SF-3  Experimental  Results 

This  section  details  the  experimental  results  using  multiple  modes  for  the  clamped 
probe  DRWG  system.  All  experimental  data  came  from  [20],  see  this  for  more  details  on 
the  experimental  setup  and  equipment  used.  In  this  experiment,  material  characterization 
measurements  were  made  at  6-18  GHz  of  a  lossy  silicon-based  magnetic  absorbing  material 
ECCOSORB  SE-3  [16]  manufactured  by  Emerson  &  Cuming  Microwave  Products,  Inc. 
Eigures  4.3  and  4.4  show  calculated  e  and  n  values  for  SE-3  respectively.  Both  figures 
show  calculations  for  1  through  5  modes,  and  also  show  experimental  results  using  a  NRW 
method  [22]. 

Eigures  4.3  and  4.4  show  that  e  and  //  calculations  using  multiple  modes  are  consistent 
with  those  using  a  single  mode  for  both  MUTs.  They  are  also  consistent  with  the  results 
using  the  NRW  method.  The  small  differences  in  results  can  be  attributed  to  measurement 
errors  and  sample  inconsistency.  While  using  higher  modes  appears  to  yield  a  better  result 
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Figure  4.3:  Calculated  e  values  of  a  clamped  DRWG  system  with  SF-3  MUT  for  1  -  5 
modes  as  well  as  experimental  results  using  a  NRW  method. 


(converges  closer  to  the  NRW  method)  for  some  frequencies  this  is  not  the  case  in  all  of  the 
frequencies.  This  leads  to  the  conclusion  that  adding  higher  order  modes  when  calculating 
6  and  yu  does  not  yield  a  significant  benefit  over  using  a  single  mode.  This  also  proves 
that  for  this  particular  system,  higher  order  modes  are  not  being  excited.  If  they  are  being 
excited,  the  amount  of  energy  coupling  is  extremely  low.  This  conclusion  is  confirmed 
by  examining  the  expected  reflection  and  transmission  coefficients  for  each  mode.  These 
expected  coefficients  (calculated  with  equation  4.28)  are  shown  in  Table  4.1,  and  show 
that  the  amount  of  energy  reflected  and  transmitted  in  the  higher  modes  is  1  to  2  orders  of 
magnitude  lower  than  that  of  the  dominant  mode. 
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|i^  and  for  SF-3  Clamped  Probe 


Figure  4.4:  Calculated  jx  values  of  a  clamped  DRWG  system  with  SF-3  MUT  for  1  -  5 
modes  as  well  as  experimental  results  using  a  NRW  method. 


Table  4.1:  Expected  reflection  and  transmission  coefficients  for  each  mode 


Mode 

IFil 

in 

1 

.6851 

.4230 

2 

.0107 

.0199 

3 

.0012 

.0005 

4 

.0054 

.0021 

5 

.0011 

.0009 
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Upon  consideration  of  the  findings  in  Figures  4.3  and  4.4,  as  well  as  Table  4.1,  it  is 
concluded  that  higher  order  modes  do  not  need  to  be  used  for  this  class  of  materials  to 
calculate  accurate  material  parameters.  Only  needing  a  single  mode  will  save  valuable 
time  and  resources,  and  will  yield  a  valuable  diagnostic  tool.  These  results  are  consistent 
to  those  found  using  the  single  probe  DRWG  system. 
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V.  Conclusion  and  Future  Work 


This  chapter  recaptures  all  work  presented,  and  details  the  results  as  well  as  presents 
possible  future  research. 

5.1  Dual  Probe  DRWG 

In  chapter  2  a  dual  probe  DRWG  system  was  analyzed.  Dual  probe  systems  have 
the  ability  to  electromagnetically  characterize  PEC  backed  materials  while  taking  both 
reflection  and  transmission  measurements  simultaneously.  Previous  DRWG  research  for 
PEC  back  material  focused  on  a  single  probe  system  which  required  a  two  thickness 
technique  to  determine  both  permittivity  and  permeability.  In  some  cases  a  two  thickness 
technique  may  not  be  possible.  This  research  removes  that  requirement  and  also  saves 
setup  time  since  a  second  setup  is  not  required  to  obtain  a  second  measurement. 

Electromagnetic  fields  inside  a  dual  probe  DRWG  system  were  presented.  Eields 
in  the  PP  region  of  the  system  were  derived,  and  the  resulting  MEIE  was  solved  using 
the  MoM.  Simulations  were  conducted  to  determine  expected  reflection  and  transmission 
coefficients  for  the  system.  The  simulations  showed  realistic  results.  Another  simulation 
was  conducted  to  determine  if  the  physical  minimum  distance  between  the  two  waveguides 
of  the  dual  probe  system  was  close  enough  to  allow  for  a  meaningful  transmission 
measurement  when  using  lossy  material.  It  concluded  that  the  transmission  would  indeed 
be  measurable  for  lossy  material.  A  final  simulation  was  conducted  using  low  loss 
material.  This  simulation  showed  that  the  transmission  coefficient  does  not  decay  as  fast 
as  it  does  in  the  case  of  lossy  material.  This  allows  for  more  freedom  when  specifying 
the  spacing  between  the  DRWG’s.  Both  simulations  for  lossy  and  low  loss  material 
exhibited  oscillation  in  the  measurements.  This  oscillation  is  attributed  to  constructive 
and  destructive  coupling  between  the  DRWG’s. 
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5.2  Multimode  Single  and  Clamped  Probe  DRWG  Systems 

In  chapters  3  and  4  electromagnetic  fields  inside  a  single  probe  DRWG  system  and 
clamped  probe  DRWG  system  were  presented  where  higher  order  evanescent  modes  were 
allowed  to  exist.  Fields  in  the  PP  region  of  the  systems  were  derived,  and  the  resulting 
MFIE’s  were  solved  using  the  MoM.  A  two  thickness  method  was  used  to  generate  two 
independent  measurements  in  the  single  probe  system.  In  order  to  solve  for  permittivity 
and  permeability,  theoretical  reflection  and  transmission  measurements  were  compared 
to  experimental  reflection  and  transmission  measurements.  An  iterative  least  squares 
technique  was  then  used  to  determine  permittivity  and  permeability. 

Results  for  permittivity  and  permeability  calculations  using  multiple  higher  order 
modes  was  presented  and  compared  to  results  using  only  the  dominant  mode.  While  results 
for  higher  order  modes  and  the  dominant  only  mode  were  consistent,  ultimately  it  was 
found  that  adding  higher  order  modes  to  the  analysis  does  not  appear  to  increase  overall 
accuracy  in  the  single  probe  or  clamped  probe  systems.  This  conclusion  was  determined 
by  comparing  these  results  to  results  gathered  for  the  same  material  using  a  NRW  method. 
Slight  differences  in  material  parameter  calculation  were  attributed  to  measurement  errors 
and  material  inconsistencies.  Finally  the  expected  reflection  and  transmission  coefficients 
for  higher  order  modes  were  used  to  confirm  that  these  higher  order  modes  are  not  being 
excited  with  enough  energy  to  influence  the  resulting  material  parameter  calculations. 

5.3  Future  Work 

This  research  has  answered  two  important  questions.  Is  there  sufficient  coupling 
between  waveguides  in  a  dual  probe  DRWG  system  to  be  used  to  characterize  lossy 
material  parameters?  Does  including  higher  order  modes  in  the  analysis  of  DRWG 
systems  increase  the  accuracy  of  material  parameter  calculations?  While  these  answers  are 
significant,  many  more  questions  in  regards  to  DRWG  systems  have  yet  to  be  answered.  A 
list  of  several  potential  topics  for  future  research  is  presented  here. 


54 


•  In  this  research,  only  simple  media  was  used  for  the  analysis  and  experimental 
measurements.  While  this  work  is  beneficial,  many  materials  used  in  electromagnetic 
systems  do  not  share  all  of  these  properties.  Many  materials  used  are  anisotropic 
where  permittivity  and  permeability  of  the  material  is  dependent  on  the  orientation 
of  the  incident  electromagnetic  field.  Attempting  to  analyze  anisotropic  material 
is  more  difficult  because  these  properties  cannot  be  simplified  to  a  single  scalar 
value.  The  resulting  permittivity  and  permeability  matrices  make  analysis  of  DRWG 
systems  more  difficult.  However,  DRWG  systems  should  be  evaluated  for  potential 
use  in  characterizing  anisotropic  materials.  It  is  recommended  that  this  research 
start  with  only  including  a  dominant  mode  reflection  and  transmission  coefficient  for 
simplicity.  However,  it  would  also  be  beneficial  to  investigate  whether  the  addition 
of  complex  media  necessitates  the  use  of  higher  order  modes  for  accurate  material 
parameter  calculation. 

•  In  this  research,  DRWG’s  were  used  based  on  their  broadband  characteristics.  Other 
waveguide  shapes  may  exist  which  could  keep  the  broadband  operation  of  the  DRWG 
yet  allow  for  higher  energy  coupling  in  the  dual  probe  system.  Such  examples 
include  a  system  in  which  a  single  DRWG  is  used,  along  with  a  compliment  shaped 
waveguide,  or  a  system  of  two  half  trough  ridged  waveguides. 

•  In  this  research,  differences  in  material  parameter  calculation  were  attributed 
to  measurement  errors  and  sample  inconsistencies.  While  these  are  reasonable 
conclusions,  it  would  be  beneficial  to  study  the  effects  of  errors  on  the  multimode 
DRWG  systems.  A  full  error  analysis  would  also  yield  higher  confidence  in  these 
calculations. 

•  In  this  research,  the  dual  probe  system  was  shown  to  be  capable  of  calculating  ma¬ 
terial  parameters  for  low  loss  media.  However,  because  reflection  and  transmission 
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coefficients  do  not  decay  very  fast  with  low  loss  media,  it  is  possible  to  have  the 
DRWG’s  farther  apart  and  still  take  meaningful  measurements.  It  was  also  noted  that 
the  optimal  placement  of  waveguides  could  change  as  operating  frequency  changes. 
Thus,  it  would  be  beneficial  to  further  study  how  the  dual  probe  system  operates  with 
low  loss  material,  as  well  as  how  error  varies  with  DRWG  spacing. 
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